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Comparison of continuous-wave C0 2 lidar calibration 
by use of Earth-surface targets in laboratory 
and airborne measurements 


Maurice A. Jarzembski and Vandana Srivastava 


Backscatter of several Earth surfaces was characterized in the laboratory as a function of incidence angle 
with a focused continuous-wave 9. 1-fxm C0 2 Doppler lidar for use as possible calibration targets. Some 
targets showed negligible angular dependence, while others showed a slight increase with decreasing 
angle. The Earth-surface signal measured over the complex Californian terrain during a 1995 NASA 
airborne mission compared well with laboratory data. Distributions of the Earth’s surface signal shows 
that the lidar efficiency can be estimated with a fair degree of accuracy, preferably with uniform 
Earth-surface targets during flight for airborne or space-based lidar. 

OCIS codes: 010.3640, 280.3640, 290.1350. 


1. Introduction 

Airborne lidar systems are being used to determine 
wind velocity 1 and to measure aerosol- or cloud- 
backscatter variability with high resolution over 
large spatial scales. 2-4 While these measurements 
are made, the lidar also often retrieves a distinct 
backscatter signal from the Earth’s surface. 1 Un- 
like atmospheric backscatter, which is highly vari- 
able both spatially and temporally, Earth-surface 
backscatter is relatively much less variant and can 
be quite predictable. Therefore routine backscat- 
ter measurements by an airborne or space-based 
lidar from designated Earth surfaces with known 
and fairly uniform backscatter properties can po- 
tentially offer lidar calibration opportunities. This 
calibration can in turn be used to obtain accurate 
atmospheric aerosol- and cloud-backscatter mea- 
surements on large spatial scales. Earth-surface 
backscatter is important, because achieving a pre- 
cise calibration factor for large pulsed lidars then 
need not rest solely on use of a standard hard target 
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(HT) procedure, which can be quite exhaustive and 
costly. Furthermore, calibration from designated 
Earth surfaces would provide an in-flight perfor- 
mance evaluation of the lidar for monitoring mis- 
alignments or laser- detector degradation. The 
potential uncertainty in a calibration factor ob- 
tained with a standard HT before field measure- 
ments may be comparable with or even higher than 
when several designated Earth surfaces with well- 
defined backscatter properties encountered during 
field measurements are used. Already, many pas- 
sive remote-sensing-satellite instruments are en- 
countering and analyzing Earth surface scattering 
at different wavelengths; e.g., the advanced very 
high resolution radiometer (AVHRR) uses ocean- 
and desert-surface scattering to calibrate its instru- 
ment occasionally. 5 Hence, with active remote 
sensing with lasers that have high-resolution data, 
calibration of a space-based lidar 6 - 7 by use of the 
Earth surfaces will be extremely useful. 

The calibration methodology with the Earth’s sur- 
face initially requires measuring backscatter of var- 
ious Earth surfaces simulated in the laboratory and 
the use of these backscatter measurements as stan- 
dards for the Earth-surface signal from airborne or 
space-based lidars. Since backscatter from the 
Earth’s surface may be retrieved at different angles of 
incidence, this work presents measured backscatter 
at various angles of incidence of different simulated 
Earth-type surfaces in the laboratory with a focused 
continuous wave (cw) C0 2 Doppler lidar at a 9.1-|nm 
wavelength. In general, Earth-surface reflectance 
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measurements have been made in the infrared, s but 
the use of lidars to characterize them and in turn use 
of the Earth's surface to calibrate lidars has not been 
made. The feasibility of this calibration methodol- 
ogy is demonstrated through a comparison of these 
laboratory measurements with actual Earth-surface 
backscatter retrieved from the same lidar during the 
NASA/Multi-center Airborne Coherent Atmospheric 
Wind Sensor (MACAWS) mission 1 on NASA’s DC8 
aircraft from 13-26 September 1995. For selected 
Earth-surface scenarios from the airborne lidar data, 
an average backscatter for the surface was estab- 
lished and the statistics of lidar efficiency with these 
selected surfaces was determined. These data were 
compared with the actual lidar efficiency determined 
with the standard calibrating hard target (CHT). 
This methodology can be carried out for other lidar 
wavelengths. 

2. Lidar Theory 

The measured signal-to-noise ratio (SNR) from a HT 
surface with backscatter p HX (sr \) (reflected power 
per steradian toward the receiver, divided by the in- 
cident power) with a cw lidar is given by the general 
expression 910 

SNR hx = T] HT if HT 7 1 p HT , (1) 

where r\ is the overall lidar system efficiency and K is 
a proportionality factor that depends on several di- 
rectly measurable lidar parameters. The transmis- 
sion efficiency T is given by T - exp( -2aL), where a 
is the atmospheric attenuation coefficient at the lidar 
wavelength over the path L to the scattering target. 1 1 
At close distances (e.g., in the laboratory), T is unity; 
however, at long distances it decreases and needs to 
be taken into account. HT can stand for a standard 
CHT or earth hard target (EHT), since both present a 
surface of scattering. Equation (1) is used with a 
CHT, such as flame-sprayed aluminum (FSA) or 
sandpaper (SND), with known p CHT to determine 
firnr- 

For a distance L from a focused cw-lidar primary 
mirror to the HT surface, K HT is given by 1012 

K ht = PirR 2 f(L ) /Bh v, (2 ) 

where range dependence f(L) is given by 

f(L) = 1/{L‘[1 + (77J? 2 /XL) 2 (1 - L/F) 2 ]}, (3) 

where hv is the laser photon energy, P is the trans- 
mitted laser power, R is the lidar-beam (1/e) 2 inten- 
sity radius at the telescope primary mirror, F is the 
distance to the center of the focal volume, and B is the 
data-system bandwidth. At the lidar-beam focal 
volume L - F, Eqs. (1) — (3) reduce to 

SNR ht = r\ UT PnR 2 Tp HT /BhvF 2 , (4) 

By use of a CHT with known p CHT , as other lidar 
parameters are measurable, t| cht can be deter- 
mined. 12 This can be done either with Eqs. ( 1)— (3) 
and measurement of SNR ( HT at various distances L 


or with Eq. (4) and measurement of SNR ( HT only at 
the focal volume. Once i] CHT is known then mea- 
sured SNR ht from any other HT surface can lead to 
Pht by inversion of either Eqs. (1M3) or Eq. (4). 
Thus Pkht f° r a variety of Earth-surface compositions 
can be characterized in the laboratory. Subse- 
quently, the characterized Earth surfaces can be used 
as calibration targets in flight to obtain r| EHT for air- 
borne or space-based lidars. Finally, once t] EHT has 
been determined, then it may be further corrected to 
^aerosol* the lidar efficiency determined with 
laboratory-generated aerosols, which results in a 
more realistic assessment of the actual lidar effi- 
ciency to be used for atmospheric measurements. 12 

3. Laboratory Experiment and Airborne Mission 

The NASA/Marshall Space Flight Center’s (MSFC’s) 
focused cw 9.1-|jLm C0 2 Doppler lidar was used for 
measuring backscatter. Details of the lidar opera- 
tion and its calibration procedure can be found else- 
where. 1 " 13 Laboratory simulation was conducted 
with a variety of homogeneous HT’s made from 
Earth-surface materials. SNR ]rr was measured 
with a spectrum analyzer (SA) as a function of inci- 
dence angle 9 (the angle between the lidar beam 
meeting the surface and normal to the surface) at the 
focal volume. Different targets were investigated to 
provide an envelope of a variety of materials that may 
be found in Earth surfaces, homogeneous (e.g., sand, 
soil, snow, water, and single-type vegetation) as well 
as heterogeneous (e.g., vegetation interspersed with 
open lard and cities). During the 1995 NASA/MA- 
CAWS mission, while atmospheric aerosol backscat- 
ter was being detected, SNR EHX was also measured 
during aircraft rolls at different ranges for various 
Earth surfaces with diverse compositions from nearly 
homogenous to quite heterogeneous mixtures. 
These measurements of SNR KHT were corrected for 
atmospheric attenuation and compared with theoret- 
ical SNF eht by use of the p EHX from simulated Earth 
surfaces measured in the laboratory. From this 
comparison, the statistics of t| eht for selected Earth 
surfaces having an average p KHX was determined and 
compared with the measured lidar efficiency by use of 
the SND CHT. 

A. Backscatter from Simulated Earth Surfaces 

Figure T shows an experimental arrangement for 
backscatter measurements of simulated Earth-type 
(a) land surfaces and (b) wet surfaces. Numerous 
materials were used to simulate land-surface HT’s: 
oak, pin ?, corn, sand, soil, and roof tar. The mate- 
rial was glued, with rubber cement, onto 8-in- 
diamete:* plexiglass disks. For the vegetation 
targets, vegetation was cut into small pieces with 
sizes of roughly several millimeters. The HT was 
attached to a rotating motor shaft to give a Doppler- 
frequency-shifted signal. It could be moved along an 
optical rail to give range response of SNR(L), or it 
could be placed at the center of the lidar-beam focal 
volume or range L = F [Fig. 1(a)]. Here 9 was 
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Earth Surface 


Fig. 1. Laboratory schematic for backscatter measurement of 
simulated Earth-type (a) land surfaces and (b) wet surfaces. Air- 
borne measurement scenario ic) of retrieval of an Earth-surface 
backscatter signal from the expanded cw-lidar beam during an 
aircraft roll. 


changed by the turning of the HT with respect to the 
lidar-beam axis. 

Simulations of wet surfaces were performed in a 
different experimental setup [Fig. 1(b)]. A rotating 
vessel containing the wet material, placed on a motor 
shaft, produced a smooth parabolic surface, owing to 
centrifugal force. The equation of the curve in two 
dimensions was geometrically determined to be 
z ~ 0.13r 2 at a constant angular speed. With the 
lidar beam directed at different spots on the parabolic 
water surface, different tp could be accurately deter- 
mined by use of the vector calculus equation tp = 
cos 1 (N • I/|N||I|), where I is the incident lidar-beam 
vector and N is the surface normal vector, which were 
geometrically measured. Other targets, such as 
melted snow simulating the Earth’s ice-capped-snow 
regions and soapy water simulating a choppy-ocean- 
water surface with white caps, were also investi- 
gated. 

With Eqs. (1) — (3), t] chx was determined from the 
measured SNR CHX (L) of two CHT’s, (SND and FSA), 
and their p CHX at 9 = 45° is given by 7.33 X 10 ^ sr 1 
and 8.28 X 10 2 sr 1 , respectively, with ~4% uncer- 
tainty. 1415 Other lidar parameter values were hv — 
2.1818 x 1 CT 20 J, P = 2.9 W ± 1%, R = 0.0305 m ± 
1.5%, F = 9.33 m ± 0.5% , and R SA = 360 kHz ± 1%, 
The uncertainty in measuring SNR in the laboratory 
with SA was —7%. The CHT’s gave t] snd — 0.17 ± 
9% and r| FSA — 0. 16 ± 9%. An average value of r| CHX 
— 0.165 ± 13% was used as the lidar system effi- 
ciency for determining the backscatter of various sim- 
ulated EHT surfaces. 

For the various land and wet Earth-surface simu- 
lations, SNR eht was measured as a function of tp at 
L = F. Then, from Eq. (4), p KH T was determined 
with the previously mentioned lidar parameter val- 
ues. Figure 2 shows p KHT as a function of tp for (a) 
land-type targets and (b) wet-type targets. Uncer- 
tainty in p EH t f° r land-type targets is — 12%, while for 
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Fig. 2. Laboratory backscatter measurements as a function of 
angle of incidence tp with the NASA/MSFC 9.1-p.m cw Doppler 
lidar from simulated Earth surfaces for (a) land-type targets and 
(b) wet-type targets. 

the wet-type targets, the uncertainty is higher. In 
general, for all the simulated EHT land-type sur- 
faces, dependence of p EH t on 9 is small. Some of the 
surfaces, such as FSA, pine needles, and uncut corn 
leaves, show an increase in p EH t f° r smaller 
whereas others, such as SND, beach sand, soil, roof 
tar, oak, and cut corn leaves, show almost no 9 de- 
pendence. These results show that there are several 
naturally occurring targets that behave like the 
standard FSA and SND CHT’s and can provide rea- 
sonable calibrating targets during flight. For the 
land-type simulated EHT surfaces, sand and soil tar- 
gets give the highest p EHX , while vegetation targets 
give the lowest. Additionally, the vegetation targets 
that were dried out had negligible change on p EHX . 
However, wetting all targets with water lowered the 
Peht dramatically for all 9. This follows the same 
trend with the wavelength-dependent reflectance 
measurement of dry and wet terrestrial surfaces in 
the near infrared, 8 thus showing that wet targets had 
lower reflectance. However, backscatter measure- 
ments of wetted targets could not be accurately quan- 
tified as a function of water content for this 
experiment, because the amount of water changed 
rapidly in time, owing to the HT angular rotation. 

For the wet surfaces [Fig. 2(b)], smooth water had 
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the lowest p EHT °f ah HT’s. This supports the pre- 
vious observation that the presence of water (p EHT ~ 
10 8 sr l ) on wetted land-type targets lowered p EHT 
as compared with dried targets ( — 10“ 3 — 10“ 1 sr -1 ). 
The presence of small ripples on the smooth water 
surface tended to increase the p EHT °f water. As the 
ripples became larger, p EHT also increased. Fur- 
thermore, the addition of soap bubbles to the smooth 
water surface greatly enhanced water p EHT by more 
than 3 orders of magnitude. High backscatter was 
also measured for snow. Despite the drastic dissim- 
ilarity between snow and soap bubbles on the water 
surface, p EHT f° r both were nearly the same. Even 
though refractive indices of water and ice are similar 
at the 9.1-p.m wavelength, backscatter appears to be 
affected solely by the surface-shape irregularities of 
snow, which shows a dramatic difference between low 
signals from smooth water surfaces and high signals 
from irregular surfaces. Also, the addition of soil to 
either soap bubbles or snow tended to further in- 
crease p EHT , which again agrees with the finding that 
the p EHT °f soils (~~10 1 sr l ) was found to be higher 
than that of snow or of water with soap bubbles 
(10 5 -10 4 sr 1 ). 

The investigated targets gave a range of possible 
variation in the p EH t that may be encountered during 
field missions from real Earth surfaces composed of 
different materials. Figure 2 shows that for the sim- 
ulated land-type Earth surfaces, backscatter at 9.1 
gm varies by a factor of —50, depending on whether 
the simulated Earth surface is purely vegetation type 
or purely nonvegetation type, or a mixture of the two. 
For mixed surfaces, if the fractions of different ma- 
terials can be estimated, then an average p E ht can he 
estimated by weighted mixing of their individual 
Pkht s - With the use of passive satellite imagery, 
such as LANDSAT, it may be possible to fractionate 
different types of surfaces encountered in the lidar- 
sampled surface area. For calibration purposes, ho- 
mogeneous surfaces would be best; but if these were 
not available, than an average p EHT of a composite 
surface could be used. 

The measurements shown in Fig. 2 were made at 
the lidar-beam focus. However, backscatter mea- 
surements made from actual Earth surfaces with an 
airborne-focused cw lidar are performed at random 
range values, depending on the aircraft altitude and 
roll. Therefore measurements of the range response 
of backscatter from the simulated Earth surfaces 
were performed in the laboratory to show that the 
measurements followed the same behavior as the the- 
ory predicted [Eqs. ( 1) — (3)] and hence to support the 
possibility that the t| kht may be determined from 
SNR EHT at any range for these targets. Figure 3 
shows the range response of HT made from beach 
sand and pine, as well as the SND CHT. The beach 
sand and pine HT’s show good agreement with lidar 
theory, which is similar to the case of SND CHT. 12 
The range response of backscatter for the other land- 
type HT’s (not shown in Fig. 3) also were in good 
agreement with the lidar theory. These results sug- 
gest that the signal from these types of Earth sur- 
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Fig. 3. Measured backscattered SNR as a function of range L for 
beach sand and pine HT's and the (SND l calibrating HT. 


faces at ranges other than L = F can be used in 
conjunction with Eqs. ( 1 )— (3) to estimate lidar effi- 
ciency. 

B. Earth -Surface Return 

Figure _(c) shows a schematic of how the NASA/ 
MSFC airborne cw-focused 9.1-fxm lidar retrieved an 
Earth-surface backscatter at the unfocused part of 
the lidar beam along with atmospheric aerosol back- 
scatter at the lidar-beam focus during an aircraft roll 
for the MACAWS mission. The lidar beam was fo- 
cused a* —54 m through a modified aircraft right- 
side-viewing germanium window in front of the 
aircraft wing. Usually, it is not overly difficult to 
distinguish backscatter between aerosols and the 
Earth’s surface. As the aircraft entered a roll, a 
small Earth-surface signal emerged at a distance L. 
As the i oil angle increased, ip decreased, which de- 
creased L; the mean signal strength increased, ac- 
cording ;o Eqs. ( 1 )— (3), and the signal could be easily 
followed for the duration of the roll. Usually this 
signal h id a Doppler-frequency shift lesser or greater 
than that of the aerosol signal, and therefore a sep- 
arate SNR for each signal could be retrieved. On 
occasion the two signals overlapped, and a separate 
SNR for both aerosol and the Earth’s surface could 
not be obtained. Generally, the Earth-surface sig- 
nal fenc ed to have a narrower spread in frequency 
than th^ aerosol signal, which has an additional 
spread, 3wing to turbulence. 

Figure 4 shows both the aerosol signal and the 
Earth-surface signal from the focused cw lidar during 
an aircraft roll, as displayed on the digital signal 
processcr in real time for 1 s of integration. The 
Earth-surface signal was from the expanded lidar 
beam and was from agricultural land located in the 
Sacramento Valley near Willows, California. With 
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Fig. 4. 1-s integrated cw-lidar spectral data taken with the digital 

signal processor of both an aerosol and an Earth-surface land 
backscatter signal at different Doppler-shift frequencies during an 
aircraft roll over Sacramento Valley near Willows, California, lo- 
cated at approximately 39 25' N. lat, 122' 15' W. long. 


propagation of a Gaussian lidar beam, the lidar-beam 
cross-sectional area near the surface was —6.4 m~ 
with a projected area of —11.6 m 2 along the surface. 
For the DC8 aircraft speed of —123 m s _1 , this cor- 
responds to Earth-surface sampling of approximately 
123-m long x —5.2 m wide, or approximately 0.15 
acres s 1 . In contrast, aerosol backscatter was 
within the lidar focal volume with a beam cross- 
sectional area of —8 X 10 5 m 2 and an integrated 
sample volume of approximately 65 m 3 . Despite 
drastic dissimilarities in beam sizes and intensities 
at the focal volume for aerosols and at the expanded 
beam at the Earth’s surface, strength in signal from 
aerosols and from the Earth’s surface can be compa- 
rable. This is due to the relatively weak aerosol 
backscatter at the high-intensity part of the lidar 
beam and the relatively strong Earth-surface back- 
scatter at the low-intensity part of the lidar beam for 
L T>> F (as shown in Fig. 6 of Ref. 12). 

The outgoing lidar beam to the Earth’s surface and 
the backscatter from the surface can undergo appre- 
ciable atmospheric attenuation. 11 The retrieved 
SNR from the Earth’s surface needs to be compen- 
sated for attenuation. The significant factors are 
given in terms of the attenuation coefficient a of ma- 
jor atmospheric components: aerosol (a = 0.024 

km 1 ), carbon dioxide (C0 2 line, a = 0.0014 km l ), 
and water (H 2 0 line, a = 0.0064 km -1 and H 2 0 
continuum, a = 0.05 km ] ). The average values for 
C0 2 and H 2 0 absorption lines were derived from line- 
by-line computations 16 while the H 2 0 continuum was 
derived from lowtran 7 with representative meteoro- 
logical conditions for the California area. 17 The av- 
erage aerosol a value was derived from the integrated 
aerosol backscatter with measured vertical profiles 
obtained from ascents and descents during the MA- 



Fig. 5. Map of region of California showing the locations where 
different Earth-surface measurements were made as illustrated by 
different symbols. 

CAWS mission. The total transmission efficiency T, 
which typically ranged from 55% to 75%, was esti- 
mated for each datum with the aircraft radar altitude 
and roll angle. 

The geographical location of the various regions 
where a backscatter signal from the Earth’s surface 
was obtained during an aircraft roll are shown Fig. 5. 
Different surfaces encountered were vegetation and 
nonvegetation agricultural areas over California’s 
San Joaquin (open inverted triangle, plus, and ant) 
and Sacramento (open square) Valleys; forest mixed 
with sparse open land over the Sierra Nevada Moun- 
tains west of Reno, Nevada (filled triangle, open di- 
amond, and open triangle), and over the Coastal 
Range Mountains in California near Monterey (filled 
circle), Santa Cruz (filled square), Salinas (filled in- 
verted triangle), and Santa Rosa (cross); and grass- 
land east of San Francisco Bay (filled diamond) and 
metropolitan areas in the South Bay area (open cir- 
cle). Although the DC8 made many rolls over the 
open ocean, no signal from a relatively calm ocean 
surface (even at L — 500 m) could be detected with 
the given output power of the lidar. This was prob- 
ably due to the water’s low backscatter (Fig. 2). 

Figure 6 shows measured SNR normalized by T 
(SNR eht = SNR/Tj at various ranges from the 
Earth-surface locations identified by the symbols 
shown in Fig. 5. Each data point represents aver- 
aging over several 1-s samples, ranging from 7 to 53 
samples. The horizontal line through each data 
point gives the L variability that is due to both roll 
angle and radar altitude, whereas the vertical line 
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Fig. (T Measured backseattered SNR normalized by transmission 
efficiency T for the various Earth surfaces shown in Fig. 5 with the 
9.1-jmm cw lidar along with comparison with lidar theory [Eqs. 
( by use of lidar parameters cited in the text for four repre- 
sentative backscatter values of various Earth-surface types. 
Variability of Earth-surface backscatter at each location is de- 
picted bv the vertical line through each data point, which gives the 
range of SNR l; u T at a given distance, while variability in L is 
depicted by the horizontal line through each data point, which 
gives the range of L that is due to variability in both radar altitude 
and roll angle. 


gives SNR EHT variability over the sampled surface. 
For some samples, SNR EHT was relatively constant, 
e.g., grass hills interspersed with some trees and 
shrubs (filled diamond) and also for metropolitan ar- 
eas (open circle) which, even though inhomogeneous, 
form a uniform mixture without significant backscat- 
ter changes between consecutive samples. Whereas 
for other regions, such as over the valley and coastal 
ranges, there is more variability of SNR KHT , caused 
by significant heterogeneity encountered as a result 
of low -backscatter vegetation interspersed with high- 
backscatter nonvegetation areas. Figure 6 also 
shows curves of predicted SNR EHT as a function of L 
with Eqs. < 1 ) — ( 3 ) for four Earth-surface types simu- 
lated in the laboratory with p p:HX = 0.08, 6 . 03 , 0.002, 
and 10 (l sr 1 (Fig. 2) corresponding to sand, soil, 
vegetation, and water with ripples, respectively. 
These curves were derived with the various lidar pa- 
rameter values during the mission as hv = 2.1818 X 
10 20 J, P = 4.4 W ± 39f, R = 0.0265 m ± 3 r /r , F = 
54.0 m ± 2 r /r , and =141 kHz it Iff during the 
MACAWS mission. The measurement error in the 
SNR on the airplane during calibration with the SND 
CHT was - 87( . The uncertainties given here for the 
mission and those given earlier for the lab measure- 
ments represent instrument-related errors. The av- 
erage lidar efficiency as measured at the beam focus 
during several calibration measurements performed 
for this two-week mission period with the SND CHT 
was Tj ( HT - 0.126 ± 189L 
Most of the data were bounded between soil and 
vegetation and follow the trend of decreasing 
SNR kht with increasing L, as predicted by theory. 


For th s metropolitan areas (open circle), which would 
be con posed mostly of high-backscatter targets such 
as coi crete-type materials, roof tar, metals, soils 
mixed with low-backscatter targets such as vegeta- 
tion, f EHT was bounded by the curves of soil and 
vegetation. The grass hills interspersed with trees 
and shrubs (filled diamond) east of the Bay area gave 
low SNR EHT , close to that predicted for vegetation. 
In this case the long grass offers a good ground cover 
over soil, so backscatter from soil was not expected. 
Forested areas interspersed with open land in the 
Sierra Nevada Mountains (filled triangle, open dia- 
mond, and open triangle) and the Coastal Range 
Mountains (cross, filled square, filled circle, and filled 
in verb ?d triangle) gave an SNR EHT between soil and 
vegetation values; however, near Salinas (filled in- 
verted triangle), more sparse vegetation on the 
eastern slopes showed the highest SNR KHT . Fur- 
thermore, SNR eht from Sacramento (open square) 
and San Joaquin (open inverted triangle, plus, and 
ant) Valleys each showed that the targets were het- 
erogeneous, mixed with vegetation and nonvegeta- 
tion types. 

The sand, soil, and vegetation targets, even though 
they were made from Alabama sources, provide a 
good envelope of possible SNR EHT from surfaces that 
may contain similar materials. To get a one-to-one 
agreement between heterogeneous surfaces sampled 
during the flight with theoretical predictions, actual 
Earth-surface samples encountered in the mission 
would be needed to determine their p EnT in the lab- 
oratory. The measured SNR KHT from a 9.1 .-pm lidar 
for different Earth surfaces encountered here ranged 
within a factor of -50, as is also seen from lab mea- 
surements in Fig. 2. 

Some comparison can be made with the Earth- 
surface signal obtained during the follow-on 1996 
MACAWS mission 1 by the calibrated 10.6-pm pulsed 
lidar. Lab characterization of p r]IT from 9.1- to 
10.6-jji n wavelengths shows similar values. 14 The 
Pkht* J t 10.6 jutm, retrieved over San Joaquin Valley 
rangec from —0.003 to 0.031 sr 1 , encompassing the 
highly variable geographical features for the valley 
and agreeing well with the 9.1-jj.m data (Fig. 6, plus 
and op?n inverted triangle) bounded by the curves of 
0.002 sr 1 < p EHT < 0.03 sr 1 . Another example 1 of 
p EHT at 10.6 fim over the Coastal Range Mountains in 
California north of San Francisco ranged from 0.0144 
to 0.0f78 sr *. This example compares well with 
the 9.]-(im data (cross) which would correspond to 
the average p EHT - 0.012 sr 1 curve. In addition, 
the p Ef [T at 10.6 pm over relatively calm ocean water 
was approximately measured to be (1.8 ± 0.2) X 
10 s * 1 . Backscatter from water surface is more 
difficult to quantify, owing to the effect of surface 
disturbances and, in addition, its wavelength depen- 
dence i 3 not well known. Even so, the 10.6-|mm value 
over oj: en ocean surface is closer to the value of — 10 8 
sr 1 for 9.1 pm for a smooth water surface rather 
than to the 10 4 sr 1 value of the perturbed water 
surfac* with bubbles. These examples show that 
Earth- surface signals can be monitored consistently 
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Fig. 7. Histograms of measured SNR normalized by the transmission efficiency 7 and f{L) from the Earth s surface at various locations 
as shown in Fig. 5: (a) Coastal Range Mountains northeast of Santa Cruz, California (solid square}; (b) agricultural land over Sacramento 

Valley near Willows, California (open square); (c) eastern slopes of the Sierra Nevada Mountains (filled triangle); (d) San Jose metropolitan 
area only at L 950 m (open circle); and (e) Coastal Range Mountains east of Salinas, California (filled inverted triangle). Histograms 
of calculated lidar efficiency r| K | 1T from the respective Earth surfaces are shown in (f)-(j). The average f° r each Earth-surface 
scenario as assessed from Fig. 6 is shown in the legend. The SND t) ( MT = 0.12 ± 1 H C A for this mission is shown within the bold vertical 
lines for comparison with t| kiit from different Earth surfaces. 


from different lidars and, therefore, may be used to 
calibrate the lidar in flight with a priori knowledge of 
the pkht- 

4. Lidar Calibration from the Earth’s Surface 

The Earth surfaces sampled in flight were, for the 
most part, heterogeneous, and the SNR EHT varied a 
fair amount. Figure 7(a)-7{e) histograms for five 
sets of data over different regions show the frequency 
of occurrence of measured SNR normalized by T and 
/XL). These histograms give the distribution of 
range-independent SNR EHT , which shows the varia- 
tion, encountered during a roll, that is due to surface 
inhomogeneity only. Figure 7(a)-7(e) shows the po- 
sition of the mode of distribution changing from 
lower- to higher-backscatter targets with the tightest 
distribution given for the metropolitan area in the 
South Bay area in Fig. 7(d), which suggests that it 
was the most uniform target with the least variability 
from each 1-s sample. 


Homogeneous or even uniform heterogeneous tar- 
gets [e.g., Fig. 7(d)] would be preferred for in-flight 
calibration, as they would give well-defined SNR EHT , 
leading to well-defined lidar efficiency. However, 
since the Earth-surface targets sampled during the 
mission were mostly heterogeneous with varying 
SNR eht , any estimation of lidar efficiency ti eht [with 
Eqs. (1) — (3)] from these data sets would also be asso- 
ciated with some variability. From the statistics of 
range-independent SNR eht shown in Fig. 7(a) — 7(e), 
the statistics of r| EHT were determined with an average 
Pkmt corresponding to each of the five regions. These 
values for p EH x of 0.007, 0.009, 0.012, 0.015, and 0.021 
sr~ 1 were estimated for the curves that best fit the 
selected data samples shown in Fig. 6. The distribu- 
tions of r| EHT are shown in Fig. 7(f) — 7(j) with the mean 
t] eht and its standard deviation for each distribution 
given in the legend. Figure 7(i), corresponding to the 
most uniformly distributed target, even though heter- 
ogeneous with the tightest distribution, gave mean 
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t) h:ht ^ 0.12 with only 36^ variation as compared with 
T| ( HT - 0.126 with 189r uncertainty with CHT SND 
used during the mission, as stated in Subsection 3.B. 
The uncertainty in the r\ KHT estimation here is domi- 
nated by the variability in the SNR h:HT from the com- 
plex heterogeneous surfaces. The most nonuniform 
Earth-surface target [Fig. 7(f)] gave mean r| EHX - 0.15 
with a variation of - 80 C A , and if an appropriately vary- 
ing pkht was known for such a target, then the esti- 
mation of t) kht would be much less uncertain. These 
distributions show that if the Earth-surface backscat- 
ter is known then the lidar efficiency can be estimated 
with a fair degree of accuracy, preferably with uniform 
Earth-surface targets during flight for airborne or 
space-based lidar. 

5. Conclusion 

This paper characterizes the lidar backscatter of sim- 
ulated Earth surfaces in the laboratory for their pos- 
sible application as calibration targets for airborne 
and space-based lidars during flights. Several dif- 
ferent simulated Earth surfaces were investigated, 
such as vegetation and nonvegetation land-type and 
wet-type surfaces. The land-type targets provide a 
fair envelope of Earth-surface backscatter for several 
materials. The backscatter of most dry Earth sur- 
faces at various angles of incidence showed only weak 
to negligible angular dependence. The SNR data 
measured in the airborne MACAWS mission over 
various complex heterogeneous surfaces lie within 
the bounds of the values obtained in the laboratory. 
Thus data from different regions could be used for 
in-flight calibration to estimate the lidar efficiency 
and the associated uncertainty. The Earth-surface- 
derived lidar efficiency obtained with the average 
Pkht °f complex heterogeneous surfaces encountered 
during the mission gave good agreement with the 
CHT-derived lidar efficiency with least error associ- 
ated with the most uniform target. Therefore, by 
measuring SNR KHT from Earth surfaces with known 
Pkht- using an airborne or space-based lidar, one can 
estimate the lidar system efficiency for atmospheric 
measurements. Obtaining absolute calibrated 
backscatter measurements with this method instead 
of signals in relative units has important advantages 
that lead to crucial information and immense re- 
search opportunities for various aspects of global 
aerosol modeling in terms of its impact on climate, 
pollution, and hydrological processes. 
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